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ABSTRACT 

Recently there has been a lot of attention focussed on a virialized halo-based approach 
to understanding the properties of the matter and galaxy power spectrum. We show 
that this model allows a natural treatment of the large and small scale redshift space 
distortions, which we develop here, which extends the pedagogical value of the approach. 



1 INTRODUCTION 



The power spectrum of the mass fluctuations in the uni- 
verse is one of the most fundamental quantities in large-scale 
structure. It is robust, but sensitive to several cosmological 
parameters such as the Hubble constant, the matter den- 
sity and of course the primordial power spectrum (usually 
parameterized by an amplitude and a slope) . While the the- 
ory behind the power spectrum in the linear regime is quite 
straightforward, analytically handling clustering in the non- 
linear regime has proven quite difficult. Recently several au- 
thors (|Ma fc Fry 200C| ; |Scljak 2000| ; [Peacock 2000| ; |Peacock 
|fc SmitF2()00| ) have developed a new way of looking at the 
non-linear power spectrum which imagi nes all the mass in 
the u niverse lies in a halo of some mass (Press & Schechter 



1974). They postulate that on large scales the halos cluster 



according to linear theory while o n small scales the power 
is dominated by the halo profiles (Neyman, Scott & Shane 
|l95£ ; Peebles 1974). Specifically most pairs of dark mat- 
i ter particles with small interparticle separations lie within 
the same halo, and thus their correlations can be predicted 
by the halo profile. While this model requires many ingredi- 
ents to be fixed by numerical experiments (typically N-body 
simulations) it provides a useful structure for thinking about 
gravitational cluste ring which g i ves insights into several out- 
standing problems ( jeljak 200C ; |Pcacock fc Smith 2000|; Sel 



ak, Burwell fc Pen 2000|; |Atrio-Barandcla fc Muckct lOOOf 



Cooray, Hu fc Miralda-Escude 2000[ ). 

The work to date has all focussed on the clustering 
of the matter or galaxies in "real space", whereas many if 
not most observations of clustering take place in "redshift 
space". It is the purpose of this work to show that red- 
shift space distortions can be handled naturally in the halo 
picture and that doing so provides insight into some well 
studied phenomena. 



2 THE HALO MODEL 

The model for n on-linear clustering is based on the Press- 
Schechter (1974) theory, in which all of the mass in the 
universe resides in a virialized halo of a certain mass. The 



calculate the power spectra and cross correlations between 
both the mass and the galaxies, given a suitable prescrip- 
tion for how galaxies populate dark matter halos. Since at 
present such prescriptions are somewhat ad hoc we shall con- 
centrate here on the mass power spectrum, though there is 
no obstacle in principle to extending the method to galaxies. 

In this model the power spectrum, P(k), is the sum of 
two pieces. The first is that due to a system of (smooth) 
halos of profile y(k) laid down with inter-halo correlations 
assumed to be a biased sampling of Pn n (k). Since the real 
space convolution is simply a Fourier space multiplication 
this contribution is 



3 (fc) = P lin (fe) 



f(u)du b{v)y{k-M) 



(1) 



where b(y) is the (linear) bias of a halo of mass M(v) and 
f(v) is the multiplicity function. The peak height v is related 
to the mass of the halo through 



(2) 



o-(M) 



where S c = 1.69 and a(M) is the rms fluctuation in the 
matter density smoothed with a top-hat filter on a scale 
R s = 3A//4?rp. Bot h b and / come from fits to N-body 
simulations. We use (sheth & Tormen 1999) 
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(3) 



(4) 



where p = 0.3 and v' = 0.707^ . The normalization constant 
A is fixed by the requirement that all of the mass lie in a 
given halo 



f(?)dv = 1 



(•>) 



We assume that the halos all have spherical profiles de- 
pending only on the mass. We neglect any substructure or 
halos-within-halos as this seems to be uni mportant for the 



power spectrum . We take the 'NFW form (Navarro, Frenk 



extensions to this theory introduced recently allow one to 



fc White 1996) 
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Figure 1. The predictions of the halo model for real and redshift 
space clustering, compared to a sequence of N-body simulations. 
The upper solid line the sum of Eqs. (|l], [lo| . The solid squares 
show the results of N-body simulations. The lower solid line is the 
prediction for the redshift space power spectrum (Eq. |l6| ) and the 
open squares the N-body results. The contribution of the 2-halo 
(dotted) and 1-halo (dashed) terms to the total is also shown. 



p(r) = 



(6) 



x(l + x) 2 

where x — r/r s is a scaled radius, thoug h our result s 
are not particularly sensitive to this choice (Seljak 2000). 
The mass of this halo is defined to be the virial mass 
M = (47r/3)5 v i r pr^ ir where we take <5 v i r = 200 and define 
the virial radius, r v i r , as the radius within which the mean 
density enclosed is 5 v i r times the background density. The 
concentration parameter c = r v ; r /r s . The profile can then 
be described by its virial mass and concentration. N-body 
sim ulati ons suggest the two are related and we follow (Seljak 
2OOC01 using 

o(M) = 10 (_) (7) 

where <r(M») = 1. 

Finally then we need y(k) which is the Fourier transform 
of Eq. (^) normalized to unit mass. The Fourier transform 
can be done analytically 



p(k) = 4irp r 3 a [ 



cosz{Ci([l + c]z) 
sinz{Si([l + c]z) 



-Ci(*)} + 
Si(z)} - 



(l + c)z 



where z = kr a and the total mass is 



M = iivpori 



log(l 



1+c 



(8) 



(9) 



Eq. (Q) is the dominant contribution to the power spec- 
trum on large scales. On small scales we are dominated by 
pairs lying within a single halo 



P 1 " 



3 W 



(27T) 



f{v)dv 



M{y) 



(10) 



We show an example of how well this formalism predicts 
the matter power spectrum at z — in Fig. |l|. We have 
chosen a particular ACDM model with Q m = 0.3, Qa = 0.7, 
h = 0.7, £l B h 2 = 0.02 and n = 1. The model has been 



normalized to t he C OBE 4-year data using the method of 
Bunn & White (1997). The agreement in the linear regime is 



good, as is to be expected. On smaller scales the model is a 
surprisingly good fit to the N-body data given the simplicity 
of the assumptions. The slight shortfall in power at higher 
k can be remedied by modifying the prescription som ewhat, 
but w e will here stick to the parameters outlined in (Seljak 
200C ) — this model is valuable more for its pedagogical value 
than as a substitute for direct calculations. 



3 REDSHIFT SPACE DISTORTIONS 

The model described above naturally lends itself to a treat- 
ment of redshift space distortions. There are two effects 
which come in when moving to redshift space. The first is 
a boost of power on large scales due to streaming of matter 
into overdense regions. The second is a reduction of power 
on small scales due to virial motions within an object. In 
the halo model the large-scale and small-scale effects can be 
separated o ut in a simple fashion. 

Kaiser (1987) first showed that on large scales one ex- 



pects an enhancement of the power spectrum in redshift 
space. In linear theory a density perturbation 8k generates 
a velocity perturbation 5 = —ikv with v parallel to k. Using 
density conservation to linear order and making the distant 
observer approximation (kr 3> 1, we shall work throughout 
in the plane-parallel limit, see (Heavens & Taylor 1995; Sza- 
lay, Matsubara & Landy 1995 ) for a large-angle formalism) 



the redshift space density contrast can be written as 



= Sr, 



1 + f» 



(11) 



where /(fi) = d log <5/d log a ~ f2 0,6 , a is the scale-factor 
and fi = f ■ k. Thus on large-scales the power spectrum is 
increased by 



dp (i + / M 2 ) 2 = i + |/+ I/ 2 



(12) 



which is approximately 1.37 in our model. 

On small scales virial motions within collapsed objects 
reduce power in redshift space. If we assume that our halos 
are isotropic, virialized and isothermal with ID velocity dis- 
persion o then the peculiar motions within the halo add a 
Gaussian noise to the redshift space radial coordinate. Once 
again the real space convolution becomes a Fourier space 
multiplication and the inferred density contrast is 



Integrating this over p gives a suppression 
V 2 y 



(13) 



(14) 



As pointed out by Peacock & Dodds ( 1994 ) however, the 
"full" effect of redshift space distortions includes both the 
enhancement and the suppression of power, and one must 
include the p dependence of both factors before doing the 
integral. Including both terms the redshift space distortion 
becomes, upon integrating over p, 
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Figure 2. The ratio of power in redshift space, compared to 
real space, from the halo model (solid line) and from N-body 
simulations (points). 
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f 2 (3 + 2y 2 )+4fy 2 ] 



(15) 



Thus to predict the redshift space power spectrum in 
the halo model we modify Eqs. (@, 0) to 

P(k) = (l + |/ + i/ 2 ) flln(fc) X 



+ 



(27T) 



/(i/)<fc- b(v)Hi(k<j)y(k; M) 



f{v)du ^^-TZ2(ka)\y(k)\ 2 



(16) 



In the first term we have broken out the enhancement due 
to halo motions, treated in linear theory, and the IZi term 
from the virial motion. Technically we should do the integral 
over v including the Gaussian from Eq. (jl3|) first and then 
integrate over jj,. We can see in Fig. ^ however that this nu- 
merically simpler approximation works quite well. Note that 
in contrast to the real-space power spectrum the clustering 
term remains significant to larger k. 

We assume that the halos are isothermal. From the mass 
within the virial radius the the ID velocity dispersion of a 
halo of mass M is 



a 2 = GM/2r viT 



G 



(^M 2 p<5 vir ) 



1/3 



(17) 
(18) 



We find that our results are almost entirely unchanged if we 
estimate a from the circular velocity interior to r s instead 

of r v ir- 

Putting the pieces together gives the lower solid line in 
Fig. [[]. Notice that the model provides an adequate descrip- 
tion of the redshift space power spectrum. The clustering 
term remains significant to smaller scales, which may ex- 
plain why perturbation theory results seem to work better 
in redshift space than in real space and why the redshift 
space power spectrum approximates the linear theory pre- 
diction over such a wide range of scales. At small-fc this 



model makes the same predictions as linear theory with a 
constant "effective" bias 



(6) = / f{v)du b[u) 



(19) 



The high-fc behavior of the model qualitatively reproduces 
the suppression seen in the N -body simulations. This is no t 
surprising since Sheth (1996) and Diaferio & Geller (1996) 
have shown that a sum of Gaussian random velocities 



weighted by the Press-Schechter (1974) mass function pro- 
vides a good description of the exponential distribution of 
velocities seen in N-body simulations. 

The model underestimates the N-body results in Fig. [I] 
in both real and redshift space. We show the prediction for 
the ratio of redshift-space to real-space power in Fig. ^ com- 
pared to the same ratio from the N-body simulations. Here 
we see that the agreement is very good over more than 2 
decades in length scale. 



4 CONCLUSIONS 

Many of the features of the power spectrum of density fluctu- 
ations in the universe can be simply understood in a model 
based virialized halos. Important ingredients in the model 
are that the halos be biased tracers of the linear power spec- 
trum and have a uniform profile with a correlation between 
the internal structure and the mass wh ich sh ould span a wide 
range sampling a a Press-Schechter (1974) like mass func- 



tion. Within this model it is easy to account accurately for 
redshift space distortions which alter the power on large and 
small scales. The effects of virialized motions within halos 
suppress the "Poisson" or 1-halo term, so the redshift space 
power spectrum traces the linear theory result to smaller 
physical scales. This effect could be responsible for the fact 
that perturbation theory is known to work better in redshift 
space than real space when compared to numerical simula- 
tions. 
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